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' ^ '■ Abstract 

The non-extensive self-consistent theory describing the thermody- 
namics of hadronic systems at high temperatures is used to derive some 
r~| ■ thermodynamical quantities, as pressure, entropy, speed of sound and 

trace-anomaly. The calculations are free of fitting parameters, and 
^ ' the results are compared to lattice QCD calculations, showing a good 

^ ■ agreement between theory and data up to temperatures around 175 

MeV. Above this temperature the effects of a singularity in the par- 
^T) ' tition function at Tq = 192 MeV results in a divergent behaviour in 

' respect with the lattice calculation. 

(N 

^ One of the challenges of nuclear physics nowadays is to understand the na- 

ture of the transition from confined to deconfined regimes of strong-interacting 
systems. In this regard the thermodynamics of hadronic medium has received 
much attention due to the possibility of describing the equation-of-state at 
temperatures just below the critical value, T^, where the system could still 
be represented in terms of hadronic degrees of freedom. Quantities as the 
trace-anomaly and the speed of sound in hadronic medium are considered 
sensitive to the phase transition. 

Hadron resonance gas (HRG) models [THH] use Hagedorn's theory [6] to 
introduce an exponentially increasing mass spectrum for hadrons. This be- 
haviour of the mass spectrum was considered to be necessary [7] to describe 
the lattice-QCD results. However, a few problems have appeared: Hage- 
dorn's theory fails to explain the experimental results from hadron-hadron 
collisions for energies above 10 GeV [8l[9]; different HRG models diverge 



1 



on the value for the Hagedorn's temperature when fitting the available in- 
formation for mass spectrum, some of them being well-above the critical 
temperature expected from lattice-QCD calculations [TO] . 

Recently the non-extensive self-consistent thermodynamical theory was 
developed predicting a limiting effective temperature, Tq, a limiting entropic 
index, g^, and a power-law mass-spectrum [TT]. This theory imposes a much 
more restrictive test to the role played by Tsallis statistics in High Energy 
Physics (HEP). Similar predictions about a limiting effective temperature 
were also proposed through a massless parton gas model [12]. 

Despite being still a rather controversial subject, the use of non-extensivity 
in HEP has received increasing interest over the last years (see for instance [H 
[^[TT | [T ^[nH2B] and references therein), and it has been shown that exper- 
imental data from high energy collisions indeed supports the conclusions 
from Refs. [Tn[T2] . In this work, another important test of non-extensivity 
is performed. Using the informations of experimental data analysis, the 
non-extensive self-consistent thermodynamics is completely described. Then 
it is used to calculate thermodynamical functions which can also be cal- 
culated by lattice QCD methods. Although power-law mass spectrum was 
discarded within extensive (Boltzmann) thermodynamics [7], in this work 
the non-extensive self-consistent predictions are investigated and compared 
to lattice-QCD calculations. It turns out that there is a good agreement be- 
tween thermodynamics and lattice calculations up to temperatures around 
175 MeV. 

The non-extensive self-consistent theory has two main ingredients [TT] : 
the non-extensive statistics and ii) the fireball self-consistency principle. This 
theory predicts the limiting temperature and the limiting entropic index 
which characterize strongly interacting systems. More relevant to the present 
work are the mass spectrum, p(m), and the partition function, Z{V,(3) de- 
rived in that theory, which are given by [TT] 

p(m) = 7m~^/^eq(/3om) (1) 

and 

\n[Z{V,(3)] = m-'p{m)e,{-f3m)dm + Z^{V,f3), (2) 

where eq{x) is the q-exponencial function 

e,(x) = [l + (g-l)x]i/(''"i), (3) 

P = 1/t, with r being the effective temperature, /3o = ^/'^o, M is the lowest 
mass for which the mass spectrum formula is compatible with the experi- 
mental information on hadronic states (see text below), and 7 is a constant. 
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The second term on the right of Eq. [2] is 



M 



m' 



^f"^ p\m)eq{— I3m)dm , 



(4) 



with p'(m) being the mass spectrum vahd below the mass M. 

Different analysis of transverse momentum distributions from high-energy 
particle collisions have shown that the prediction of the limiting tempera- 
ture and of the limiting entropic index are in accordance with experimental 
data [HI El [IMS], resulting in =(60.7±0.5) MeV, Qo =1.138±0.006 and 
7 = (5 ± 3) X 10"'^ Gev'^/^. Also, it was shown that and Qo are related 
through the equation 



where c is a constant depending on the energy transfer between the source 
and its surroundings [TB]. Since r — )• Tq as g ^ 1, Tq is considered to be the 
Hagedorn's temperature. From a systematic analysis of experimental data it 
was found that To = Th =(192±15) MeV and c =-(950±10) MeV [I3]. Ob- 
serve that the Hagedorn's temperature obtained in this way is in agreement 
with the lattice-QCD resuhs p71l28]. 

The mass spectrum given by Eq. [T]has been compared with the observed 
data. It was shown that it gives an excelent description of the available data 
for masses going from the pion mass, m^r, up to m = 2.5 GeV with values for 
To and Qo similar to those obtained in the analysis of transverse momentum 
distributions [15]. This result allows one to identify M = m^r, therefore it 
is reasonable to conclude that p'(m) = 0, and consequently Zi{V, (3) = 0. 
Hence the non-extensive self-consistent theory can completely describe the 
experimental information about the hadronic system formed in high-energy 
collisions and about the observed hadron-mass spectrum. 

In the following Greek letters are used to designate thermodynamical 
quantities which are calculated according to the non-extensive self-consistent 
theory and Latin letters to designate those quantities calculated by lattice- 
QCD. From equations [Hand [2] it is possible to calculate all thermodynamical 
functions. The fact that [{qo - - /3o)M]-^/(^~^) > 1 in the range of 
temperatures of interest is used to simplify the expressions for the thermo- 
dynamical functions presented below. It is useful to write Eq. |2] in terms 
of the Gauss' hypergeometric function, 2-^1, resulting (using ( to indicate 
partition function) [TT] 




(5) 
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where 

" 1 1 g -1 



g-l'g-l'g-1' (g - - /3o)M 

(7) 

From the equation above it is possible to get the pressure, 11, through the 
relation 

n = ^ln[C(\/,/3)] (8) 

resulting in 

n = /3-'/S^i(g,/3). (9) 

Using the above expression for 11 and the relation 

dli , , 

(10) 

the entropy density cr can be obtained, resulting 

cy=\fi-''l\7M.fi)- (11) 
The energy density, e, can be calculated from Eq. [6] by 



resulting in 



e = ^^ln[C(\/,/3)], (12) 



^ = ^r'/'2-^i(g,/3). (13) 



Other interesting quantities to be calculated are the so-called trace-anomaly, 
a(T), which is considered as an indication of phase transition, and the specific 
heat at constant volume. The trace-anomaly is calculated by 

e-m _ (9 / n 



-(T)^ — =r-y-y (14) 

and the specific heat at constant volume by 

(15) 

Also, from the entropy density and specific heat it is possible to calculate the 
speed of sound in the hadronic medium, z/, by 

.2 ^ 



v' = -. (16) 
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Figure 1: Non-extensive functions: normalized entropy density (a) and spe- 
cific heat (b) as a function of the effective temperature. 



The results for the normalized entropy and specific heat are plotted in 
Fig. [H It is possible to observe that both quantities present a fast increase 
as the temperature approaches the critical value, To, which is related to the 
existence of a singularity in the partition function [TT]. This singularity is 
the cause of the limiting effective temperature in the thermodynamical non- 
extensive self-consistent theory. 
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Figure 2: The function In Z as a function of k. The fiat region between 
k =0.1 and k = 0.4 shows that InZ is proportional to k^. 

In order to compare the non-extensive calculations with the lattice calcu- 
lations it is necessary first to find a map between the effective temperature, 
r, and the lattice-QCD temperature, T. Equation |5] establishes a connection 
between the critical temperature, Tq, and the critical effective temperature. 
To- This relation, however, cannot be used for arbitrary temperatures. The 
map between T and r is expected to be a function r(T) with the following 
properties: 





5 



1. r(To) = To] 

2. r(0) = 0; 

3. r(Ti+T2) =r(Ti) + r(r2). 

The third property will be discussed below. 

The function satisfying all those conditions simultaneously is 

T{T) = kT^^T, (17) 
where the constant k = To /To has been defined. 




T(MeV) 
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Figure 3: Normalized pressure (a) and normalized energy density (b) as 
a function of the temperature. Thick line represents the theoretical results. 
Thin lines represent the upper and lower limits of the theoretical calculations 
considering the uncertainties in the parameter 7 of the mass spectrum. Full 
circles represent the results of lattice QCD from Ref. [22]. Empty symbols 
are the results of lattice calculation from Ref. [30], with circles indicating p4, 
N£=4, squares indicating p4, Ni=6. 



Note that Eq. [T7] determines also relations between other thermodynam- 
ical functions. In fact the partition function is transformed in such a way 
that 

\nZ=k-^\nC (18) 

for k around the value To/To, as can be observed in Fig. [21 

With these relations established, it is now possible to compare the results 

from the non-extensive self-consistent theory with lattice calculations. 

In the present work three sets of lattice data are used for comparison 

with the theoretical results, namely, that of S. Borsanyi et al. [29] (LI), 
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and those from M. Cheng et al [30] (L2). The former is a 2+1 staggered 
quarks calculation with physical quark masses, stout action and A^^ in the 
continuum limit, while the calculations in the latter set are obtained with 2+1 
staggered quarks with almost physical masses and p4 action for Nt =4 and 
Nt =6 . These three sets are not intended to provide a complete description 
of lattice-QCD results, but are representative of the state-of-art in the field. 

The results for the normalized pressure are shown in Fig. EJ^a), and are 
compared to LI results. Considering the uncertainties in the parameter ^ it 
is possible to observe a fair agreement between theory and lattice calculations 
up to T ~185 MeV, the deviation thereof being due to the singularity of the 
partition function at the critical temperature. 

In Fig. El^b) the normalized energy density is compared to L2 data. In 
this case there is a larger deviation of lattice results in respect with the 
theory, although the overall shape of curve and data are quite similar in 
the region up to T =200 MeV. This deviation may be attributed, at least 
partially, to the use of non-physical quark masses. In fact, it is known that the 
critical temperature from lattice calculations decreases as the quark masses 
approach the physical values [3|[HT] . It was shown [5] that the use of physical 
quark masses would result in a shift in the critical temperature towards low 
temperatures of about 20 MeV. It is possible to see in Figs. [3] and H] that 
such a shift would lead to a good agreement between L2 data and theory. 

In Fig. m^a) the trace-anomaly is reported. One can see that up to 
T ~175 MeV the theoretical calculations fall between the two kinds of 
lattice-QCD results, being systematically above the results from L2 group 
and systematically below the results from LI. Considering the uncertainties 
in the parameter 7, it is possible to find an agreement with the calculations 
from Ref. [2H], but it is not possible to discard the agreement with those 
from Ref. [5U] . Despite the same critical temperature, the deviation between 
thermodynamical and lattice calculations starts at lower temperatures in the 
case of the trace-anomaly, reflecting the higher sensitivity of the latter to the 
critical temperature. 

The speed of sound in hadronic medium is plotted in Fig. Ht^b) and com- 
pared with LI calculations. Now a nice agreement is found between lattice 
calculations and theoretical results due not only to the larger uncertainties in 
lattice data but also to the fact that the speed of sound is independent of 7, 
which is the most uncertain parameter in the theory. As it was expected [32] 
the speed of sound goes to zero at the critical temperature. 

^One must keep in mind that the parameters in the non-extensive self-consistent the- 
ory are not free of correlations, therefore the uncertainty range presented here may be 
somewhat overestimated. 
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In it was shown that finite-size effects can give rise to non- 

extensivity and that the ratio r/T should be proportional to exp{—S/C), 
where S is the entropy and C the heat capacity. Using Cy for C and ob- 
serving the practically constant speedy of sound in Fig. ID^b), we see that 
the linear property for r(T) used here is in qualitative agreement with the 
results in [331131] • 

It is not the objective of this work to determine which of the lattice-QCD 
calculations are in better agreement with the non-extensive self-consistent 
theory. For this purpose more precise determination of the parameters and 
more accurate calculations are required. Such an accuracy may be achieved, 
for instance, by considering distinct mass spectra for bosons and fermions. 
However, the results presented here show that there is a better agreement 
between the theoretical results and LI or L2 than between LI and L2. 



a 




(a) (b) 

Figure 4: Trace-anomaly (a) and squared sound speed (b) as a function of 
the temperature. Symbols and lines have the same meaning as in Fig. [3] 

The results obtained so far indicate that the non-extensive self-consistent 
theory is able to directly connect experimental information to lattice-QCD 
data since it describes on one hand the experimental results for pr-distributions 
and the observed hadronic states, and on the other hand the thermodynami- 
cal functions as calculated in lattice-QCD. In this sense the theory represents 
an interesting opportunity to understand the phase transition from confined 
to deconfined regimes. 

It is worth to mention that HRG models can describe quite well the 
lattice data [35]. The main difference between the approach presented here 
and HRG results is that the non-extensive self-consistent theory is directly 
based on experimental data for p^^-distributions, while the latter approach 
lacks this base. 

Concluding, in this work several thermodynamical functions of interest 
in the study of excited hadronic systems are calculated according to the 
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non-extensive self-consistent thermodynamical theory proposed in Ref. 
A connection between the non-extensive functions and the corresponding 
lattice-QCD quantities is estabhshed. This allows the comparison of the re- 
sults obtained by using the non-extensive thermodynamics to those obtained 
in lattice-QCD calculations. 

The parameters T, q and 7 used in the thermodynamical theory are those 
determined directly from experimental p^-distributions or from hadron spec- 
trum [15], therefore no free parameters are used to fit the results of lattice 
QCD. The comparison shows the agreement between theory and lattice cal- 
culations for all thermodynamical functions analyzed here. Also, it shows 
that a power-law mass spectrum is consistent with lattice results within the 
non-extensive thermodynamics. 

With these results, the non-extensive self-consistent theory shows to be 
an important tool for the study of thermodynamical properties of hadronic 
medium. In fact, it completely describes the experimental information from 
high-energy collisions and from the observed mass spectrum, and on the 
other hand it is in agreement with the lattice-QCD calculations for the bulk 
thermodynamic functions. 
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